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Abstract 

In this article we show that the three-particle GHZ theorem can be 
reformulated in terms of inequalities, allowing imperfect correlations due 
to detector inemcencies. We show quantitatively that taking into account 
these inefficiencies, the published results of the Innsbruck experiment sup- 
port the nonexistence of local hidden variables that explain the experi- 
mental results. 

The issue of the completeness of quantum mechanics has been a subject of in- 
tense research for almost a century. Recently, Greenberger, Home and Zeilinger 
(GHZ) proposed a new test for quantum mechanics based on correlations be- 
tween more than two particles pj]. What makes the GHZ proposal distinct from 
Bell's inequalities is that they use perfect correlations that result in mathemat- 
ical contradictions. The argument, as stated by Mermin in goes as follows. 
We start with a three-particle entangled state 

W = ^(l+>l|+>2|-)3 + |-)l|-> 2 |+>3). 

This state is an eigenstate of the following spin operators: 

A = <JlxO'2yd'3yi B = <7ly<72x03s) 

C = <J\y<72y&3x, D = b\ x OlxGZx- 
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From the above we have that the expected correlations E(A) = E(B) = 25(C) = 
1. However, D = ABC, and we also obtain that, according to quantum me- 
chanics, E(D) = 25(ABC) = —1. It is easy to show that these correlations yield 
a contradiction if we assume that spin exist independent of the measurement 
process. 

GHZ's proposed experiment, however, has a major problem. How can one 
verify experimentally predictions based on perfect correlations? This was also 
a problem in Bell's original paper. To "avoid Bell's experimentally unrealistic 
restrictions", Clauser, Home, Shimony and Holt ^ derived a new set of inequal- 
ities that would take into account imperfections in the measurement process. A 
main purpose of this article is to derive a set of inequalities for the experimen- 
tally realizable GHZ correlations. We show that the following four inequalities 
are both necessary and sufficient for the existence of a local hidden variable, or, 
equivalently 0, a joint probability distribution of A, B, C, and ABC, where 
A,B, C are three ±1 random variables. 

- 2 < E(A) + 25(B) + E(C) - E(ABC) < 2, (1) 

- 2 < -E(A) + E(B) + 25(C) + E(ABC) < 2, (2) 

- 2 < E(A) - 25(B) + 25(C) + E(ABC) < 2, (3) 

- 2 < E(A) + 25(B) - 25(C) + E(ABC) < 2. (4) 

For the necessity argument we assume there is a joint probability distribution 
consisting of the eight atoms abc, . . . , abc, where we use a notation where a is 
A = l,aisA = -l, and so on. Then, E(A) = P(a) - P(a), where P(a) = 
P(abc)+P(abc)+P(abc)+P(abc), and P(a) = P(abc)+P(abc)+P(abc)+P(abc), 
and similar equations hold for E(B) and E(C). Next we do a similar analysis of 
25( ABC) in terms of the eight atoms. Corresponding to (|l|), we now sum over 
the probability expressions for the expectations F — E(A) + E(B) + E(C) — 
E(ABC), and obtain 

F = 2[P{abc) + P(abc) + P{abc) + P{abc)} 

-2[P(abc) + P(abc) + P(abc) + P{abc)}. 

Since all the probabilities are nonnegative and sum to < 1, we infer (0) at once. 
The derivation of the other three inequalities is similar. To prove the converse, 
i.e., that these inequalities imply the existence of a joint probability distribution, 
is slightly more complicated. We restrict ourselves to the symmetric case P(a) = 
P(b) = P{c) = p, P(ABC = l) = q and thus 25(A) = 25(B) = 25(C) = 2p - 1, 
25(ABC) = 2q — 1. In this case, (0) can be written as < 3p — q < 2, while the 
other three inequalities yield just < p + q < 2. Let x = P(abc) = P(abc) — 
P(abc), y = P(abc) = P(abc) = P(abc), z = P(abc) and w = P(abc). It is 
easy to show that on the boundary 3p — q defined by the inequalities the values 
x = 0, y = |,2 = 0, w = 1 — q define a possible joint probability distribution, 
since 3x + 3y + z + w = 1. On the other boundary, 3p = q + 2 a possible joint 
distribution is x = , y — 0, z — q, w — 0. Then, for any values of q and p 
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within the boundaries of the inequality we can take a linear combination of these 
distributions with weights 3p ^~ g and 1 — 3p ^~ g and obtain the joint probability 
distribution, x = (1 - = ^f, z = (1 - ^2)5, w = ^(1 - g), 

which proves that if the inequalities are satisfied a joint probability distribution 
exists, and therefore a local hidden variable as well. The generalization to the 
asymmetric case is tedious but straightforward. 

The correlations present in the GHZ state are so strong that even if we 
allow for experimental errors, the non-existence of a joint distribution can still 
be verified. Let 

(i) E(A) = E(B) = E(C) > 1 - e, 

(ii) £(ABC) < -1 + e, 

where e represents a decrease of the observed correlations due to experimental 
errors. To see this, let us compute the value of F defined above, F — 3(1 — 
e) — (— 1 + e). But the observed correlations are only compatible with a local 
hidden variable theory if F < 2, hence e < \. Then, in the symmetric case, 
there cannot exist a joint probability distribution of A, B and C satisfying (i) 
and (ii) if e < 1/2. 

We will give an analysis of what happens to the correlations when the de- 
tectors have efficiency d € [0, 1] and a probability 7 of detecting a dark photon 
within the window of observation when no real photon is detected. Our analysis 
will be based on the experiment of Bouwmeester et al. ||. In their experiment, 
an ultraviolet pulse hits a nonlinear crystal, and pairs of correlated photons are 
created. There is also a small probability that two pairs are created within a 
window of observation, making them indistinguishable. When this happens, by 
restricting to states where only one photon is found on each output channel to 
the detectors, we obtain the following state, 

-^I+)t(|+)i|+) 2 |-)3 + |-)i|-) 2 |+) 3 ), 

where the subscripts refer to the detectors and + and — to the linear polarization 
of the photon. Hence, if a photon is detected at the trigger T (located after a 
polarizing beam splitter) the three-photon state at detectors D\ , D2 , and D 3 is 
a GHZ-correlated state (see FIG. 1). 

We will assume that double pairs created have the expected GHZ correla- 
tion, and the probability negligible of having triple pair produtions or of having 
fourfold coincidence registered when no photon is generated. (Our analysis is 
different from that of Zukowski who considered only ideal detectors.) Two 
possibilities are left: i) a pair of photons is created at the parametric down 
converter; ii) two pairs of photons are created. We will denote by P1P2 the 
pair creation, and by pi-..p4 the two-pair creation. We will assume that the 
probabilities add to one, i.e. P {pi. . .pi) + P (P1P2) = 1- 

We start with two photons. P1P2 can reach any of the following combinations 
of detectors: TD U TD 2l TL> 3 , D X D U DiD 2 , D X D 3 , D 2 D 2 , D 2 D 3 , D 3 D 3 , TT . 
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FIG. 1 

Figure 1: Scheme for the Innsbruck GHZ experiment. The GHZ correlations 
are obtained when all detectors T,Di,D2, and D3 register a photon within the 
same window of time. 

For an event to be counted as being a GHZ state, all four detectors must fire 
(this conditionalization is equivalent to the enhancement hypothesis). We take 
as our set of random variables T,Di,D 2 ,D 3 which take values 1 (if they fire) 
or (if they don't fire). We will use t,di,d,2,ds (i, g?i, c^, eta) to represent the 
value 1 (0). We want to compute P (tdid 2 d 3 | P1P2) , the probability that all 
detectors T, Di, D2, D3 fire simultaneously given that only a pair of photons 
has been created at the crystal. We start with the case when the two photons 
arrive at detectors T and D3. Since the efficiency of the detectors is d, the 
probability that both detectors detect the photons is d 2 , the probability that 
only one detects is 2d(l — d) and the probability that none of them detect is 
(1 — d) 2 . Taking 7 into account, then the probability that all four detectors fire 
is 

P {tdxdkdi I V1P2 = TD 3 ) =~f 2 (d + 7(1 - d)) 2 , 

where P1P2 = TD3 represents the simultaneous (i.e. within a measurement 
window) arrival of the photons a the trigger T and at D3. Similar computations 
can be carried out for pip2 = TD\, TD2, D1D3, D1D2, D2D3. For p±p2 = DiDi 
the computation of P (idi^cfa | P1P2 — DiDi) is different. The probability that 
exactly one of the photons is detected at Di is d(l — d) and the probability that 
none of them are detected is (1 — d) 2 . Then, it is clear that 

P (tdnkda I P1P2 = DiDi) = d (1 - d) 7 3 + (1 - d) V, 

and we have at once that 

P{td 1 d 2 dz\p 1 p 2 ) = 6 7 2 (d + 7 (l-d)) 2 

+4 7 3 (l-rf) (d + 7). 
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We note that the events involving P (td\d2d 3 \ P1P2) have no spin correlation, 
contrary to GHZ events. 

We now turn to the case when four photons are created. The probability 
that all four are detected is d 4 , that three are detected is 4d 3 (l — d), that two 
are detected is 6d 2 (l — d) 2 , that one is detected is 4d(l — d) 3 , and that none is 
detected is (1 — d) 4 . If all four are detected, we have a true GHZ-correlated state 
detected. However, one can again have four detections due to dark counts. We 
will write Pi-.-Pa — GHZ to represent having the four GHZ photons detected, 
and Pi-.-Pa = GHZ as having the four detections as a non-GHZ state. We can 
write that 

P (tdid-zda I Pi-Pa = GHZ) = d 4 + 7 (1 - d) d 3 (5) 

and 

Pfaxdtdz I p 1 ...pA = GHZ) = 3jd 3 (l-d)+6~f 2 d 2 (l-d) 2 +4j 3 d(l-d) 3 +-f 4 {l-d) 4 

The last term in (||) comes from the unique role of the trigger T, that needs to 
detect a photon but not necessarily one that has a GHZ correlation. 

How do the non-GHZ detections change the GHZ expectations? What is 
measured in the laboratory is the conditional correlation P(SiS2S3 | tdid 2 d 3 ), 
where Si, S2 and S3 are random variables with values ±1, representing the spin 
measurement at D\,D% and D 3 respectively. We can write it as 

^(SxS^ I td 1 d 2 d 3 k GHZ)P{GHZ) 

hj 1313203 ta\(l 2 a 3 ) = ^=^= • 

V 1 ' P(GHZ) + P(GHZ) 

since for non-GHZ states we expect a correlation zero for the term 

E (S1S2S3 I td!d 2 d 3 hGHZ) P(GHZ) 
P(GHZ) + P(GHZ) ' 

Neglecting terms of higher order than -f 2 , using 7 <C d, and P(pip2) ^> P(pi—Pa), 
we obtain, from P{GHZ) = 6P( Pl p 2 h 2 d 2 +3P(pi...pAh(l-d)d 3 and P (GHZ) = 
P(pi...pA) [d 4 + 7 (1 - d) d 3 ] , that 

p ^ st! ,,, , ^ j?(S 1 S 2 S3 I tdMkGHZ) 

E (S 1 S 2 S 3 I td,d 2 d 3 ) = =- ■ - — — . (6) 

+ D P(pi...P4) d? _ 

This value is the corrected expression for the conditional correlations if we 
have detector efficiency taken into account. The product of the random vari- 
ables S1S2S3 can take only values +1 or — 1. Then, if their expectation is 
E (S1S2S3 I td\d 2 d 3 ) we have 

p( Sl s 2 S3 = 1 1 tdxdxk) = i + ^SiS^ I^^W 

The variance a 2 for a random variable that assumes only 1 or —1 values is 
4P(1) (1 — P(l)) . Hence, in our case we have as a variance 

g 2 = 1- [P(SiS 2 S 3 \td 1 d 2 d 3 )f . 
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Figure 2: Contour plot of the correlation as a function of 7 and d. The region 
where the correlation is 0.92 defines a region for the parameters 7 and d that is 
compatible with the Innsbruck results. 

We will estimate the values of 7 and d to see how much E (S1S2S3 | tdidzdz) 
would change due to experimental errors. For that purpose, we will use typical 
rates of detectors J?j] for the frequency used at the Innsbruck experiment, as well 
as their reported data [^). First, modern detectors usually have d = 0.5 for the 
wavelengths used at Innsbruck. We assume a dark-count rate of about 3 x 10 2 
counts/s. With atime window of coincidence measurement of 2x 10~ 9 s, we then 
have that the probability of a dark count in this window is 7 = 6 x 10~ 7 . From 
[|| we use that the ratio P{p\P2)/ P{pi---P2) is on the order of 10 10 . Substituting 
this three numerical values in (||) we have i?(SiS2S3 | tdid^dz) = 0.9. From 
this expression it is clear that the change in correlation imposed by the dark- 
count rates is significant for the given parameters. However, it is also clear that 
the value of the correlation is quite sensitive to changes in the values of both 
7 and d. We can now compare the values we obtained with the ones observed 
by Bouwmeester et al. for GHZ and GHZ st ates || , In their case, they claim 
to have obtained a ratio of 1 : 12 between GHZ and GHZ states. In this 
case the correlations are E (S1S2S3 | idicfocfa) — 0.92. It is clear that a detailed 
analysis of the parameters would be necessary to fit the experimental result to 
the predicted correlations that take the inefficiencies into account, but at this 
point one can see that values close to an experimentally measured 0.92 can be 
obtained with appropriate choices of the parameters d and 7 (see FIG. 2). This 
expected correlation also satisfies 

£(SiS 2 S 3 I tdtdzda) > 1 - 5. (?) 
This result is enough to prove the nonexistence of a joint probability distribu- 
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tion. We should note that the standard deviation in this case is 



(7 = V(l + 0.92) (1 - 0.92) = 0.39. (8) 

As a consequence, since 0.92 — 0.39 = 0.53, the result 0.92 is bounded away from 
the classical limit 0.5 by more than one standard deviation (see FIG. 3). 

We showed that the GHZ theorem can be reformulated in a probabilistic way 
to include experimental inefficiencies. The set of four inequalities (0)-(Q) sets 
lower bounds for the correlations that would prove the nonexistence of a local 
hidden- variable theory. Not surprisingly, detector inefficiencies and dark-count 
rates can change considerably the correlations. How do these results relate to 
previous ones obtained in the large literature of detector inefficiencies in exper- 
imental tests of local hidden-variable theories. We start with Mermin's paper 
[[}), where an inequality for F similar to ours but for the case of n-correlated 
particles is derived. Mermin does not derive a minimum correlation for GHZ's 
original setup that would imply the non-existence of a hidden- variable theory, 
as his main interest was to show that the quantum mechanical results diverge 
exponentially from a local hidden-variable theory if the number of entangled 
particles increase. Braunstein and Mann || take Mermin's results and estimate 
possible experimental errors that were not considered here. They conclude that 
for a given efficiency of detectors the noise grows slower than the strong quantum 
mechanical correlations. Reid and Munru obtained an inequality similar to 
our first one, but there are sets of expectations that satisfy their inequality and 
still do not have a joint probability distribution. In fact, as we mentioned ear- 
lier, our complete set of inequalities is a necessary and sufficient condition to 
have a joint probability distribution. 

We have used an enhancement hypothesis, namely, that we only counted 
events with all four simultaneous detections, and showed that with the coinci- 
dence constraint a joint probability did not exist in the Innsbruck experiment. 
Enhancement hypotheses have to be used when detector efficiencies are low, 
but they may lead to loopholes in the arguments about the nonexistence of 
local hidden-variable theories. Loophole-free requirements for detector ineffi- 
ciencies are based on the analysis of Jll]] for the Bell case and for Jl^] for the 
GHZ experiment without enhancement. However, in the Innsbruck setup en- 
hancement is necessary, as the ratio of pair to two-pair production is of the 
order of 10 10 §. Until experimental methods are found to eliminate the use of 
enhancement in GHZ experiments, no loophole-free results seem possible. 

FIG. 3 shows the number of standard deviations, as computed above, by 
which the existence of a joint distribution is violated. We can see that if we 
change the experiment such that we reduce the dark-count rate to 50 per sec- 
ond, instead of the assumed 300, a large improvement in the experimental re- 
sult would be expected. Detectors with this dark-count rate and the assumed 
efficiency are available 0. We emphasize that there are other possible exper- 
imental manipulations that would increase the observed correlation, e.g. the 
ratio P(pip2) / P(pi---P2), but we cannot enter into such details here. The point 
to hold in mind is that FIG. 3 provides an analysis that can absorb any such 
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Figure 3: Number of er's separating any observed correlation and the critical 
boundary 0.5. The square represents the reported correlation for the Innsbruck 
experiment, and the diamond represents the expected correlation if the dark 
count is reduced to 50 counts/s. 

changes or other sources of error, not just the dark-count rate, to give a measure 
of reliability. 

We would like to thank Prof. Sebastiao J. N. de Padua for comments, as 
well as the anonymous referees. 
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